The di-Q in the extended quark derealization, color screening 

model 



Hourong Pang 

Department of Physics, Nanjing University, Nanjing, 210093, P. R. China 

Jialun Ping 

CN ! Department of Physics, Nanjing Normal University, Nanjing, 210097, P.R. China; 

O 

Center for Theoretical Physics, Nanjing University, Nanjing, 210093, P.R. China 

C , 

^ , Fan Wang 

i— » . 

CNJ ; Center for Theoretical Physics and Department of Physics, 

Nanjing University, Nanjing, 210093, P. R. China 



X 



T. Goldman 



> 
OO 

Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545, USA 

O ' 

g ; Abstract 

The fifi (SIJ=-6,0,0) dibaryon state is studied with the extended quark-delocalization color- 
screening model including a it meson exchange tail, which reproduces the properties of the deuteron 
quantitatively. We find the mass of the di-O, to be about 45 MeV lower than the $7 — 17 threshold. 



The effect of channel coupling due to the tensor force has been calculated and found to be small 
in this case. We have also studied the effect of other pseudoscalar meson exchanges and sensitivity 
to the short-range cutoff radius, ro, for the meson exchanges. 
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I. INTRODUCTION 



Quantum ChromoDynamics (QCD) is believed to be the fundamental theory of the strong 
interaction. However it is difficult to use to calculate the low energy properties of complicated 
quark-gluon systems directly, such as the NN interaction and the structure of multi-quark- 
gluon systems. Various quark models afford a framework to understand the physics in these 
cases. 

In the last 30 years, various QCD inspired models have been developed with both successes 
and failures in explaining low energy hadron physics. To mention a few but certainly not 
all, there are the MIT bag model 0, cloudy bag model 0, Friedberg-Lee nontopological 
soliton model|3|], Skyrme topological soliton model||], the constituent quark modelf5|, etc. 
Different models use quite different effective degrees of freedom, which might be indicative 
of the nature of low-energy QCD. Under different approximations one can "derive" these 
models from QCD|J. Recently there has been a hot debate regarding the proper effective 
degrees of freedom of the constituent quark model [0]. 

We performed a phenomenological study of the baryon interactions with three con- 
stituent quark models, the Glozman-Riska-Brown quark-meson coupling model|8[], the Fuji- 
wara quark-gluon-meson coupling model - one version of the Manohar-Georgi chiral quark 
model||, and the quark derealization, color screening model(QDCSM) developed by our- 
selves - a modified version of the de Rujula-Georgi-Glashow-Isgur quark gluon coupling 
model f|T0], |TT| , [T2|| , and found that, in about 2/3 of the channels formed from the octet and 
decuplet baryons, these three models gave qualitatively the same effective baryon-baryon 
(B-B) interactions [[yj. These constituent quark models have the same long range Goldstone 
boson exchange interaction, (except for the QDCSM,) but have different mechanisms for 
the intermediate and short range interactions when they are used to study B-B interactions. 
Present data on baryon spectroscopy and B-B interactions seem to be insufficient to distin- 
guish between these models and they all seem to have an approximate QCD basisQ. A new 
generation of hadron spectroscopy should be helpful in disentangling these hadronic models. 

There are good reasons, both theoretically and phenomenologically, to believe that a two 
body confinement potential might be a good approximation for single hadrons. However 
there is no compelling reason to believe it is also a good approximation for multihadron 
systems. Dibaryons provide a good testing ground for different confinement mechanisms 
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and baryon interaction models. Being the unique, stable, B=2 system, and its size, binding 
energy and quadrupole moment having been measured very precisely, the deuteron plays a 
vital role in the development of NN interaction models. Since the H-particle, a dibaryon 
with strangeness S = —2, isospin 1=0 and spin J=0, was first predicted by the MIT bag 
model in 1977|fL4{|, tremendous efforts have been made both experimentally and theoretically 
to search for it. Up to now there is no experimental evidence for the H. 

Other dibaryon candidates have also been discussed in the literature. An S=0, J p = 0~, 
M~2.06 GeV, T~0.5 MeV, 1=0 or 2 dibaryon, df, was a hot topic in the 1990's||, [16 



The S=0, 1=0, J p = 3 + dibaryon, d*, has also been a long standing topic; almost all quark 
models predict there should be attraction in this channel but with differing strengths in 



the different models^ 0> EDI 0> 13- In 1987 > we showed that the S=-3, 1=1/2, J=2 
dibaryon state might be a narrow resonance in a relativistic quark model [21]. In 1990, 
V. B. Kopeliovich et a/.|22| predicted that there are strong interaction stable dibaryons 
with high strangeness, such as an S=-6, di-Q within the flavor SU(3) Skyrmion model. A 
general survey of dibaryon states carried out with both the nonrelativistic QDCSM and the 
relativistic version, using an adiabatic calculation, found few high strangeness states with 
masses around the lowest threshold [Tj| |19|]. Zhang et al. ||23|| predicted that the Qfl 1=0, 
J=0 dibaryon is strong interaction stable using a chiral SU(3) quark model. 

In this study, we use the extended QDCSM |T2|] to calculate the eigen energy of the di-f2 
state SI J = —6, 0, 0. This model allows the quark system to choose its most favorable 
configuration through its own dynamics in a larger Hilbert space than the other models, 
and takes into account the possible difference of the confinement interaction inside a single 
baryon and between two color singlet baryons. In particular, this is the unique model that 
gives an explanation of the similarity between molecular and nuclear forces. It also explains 
why nuclei are accurately viewed as a collection of nucleons rather than a single big bag 
with 3 A quarks. 

Here, we recalculate the QDCSM result because the extended QDCSM reproduces the 
deuteron properties well quantitatively, unlike the QDCSM, and because preliminary NN 
phase shift calculations using it fit the NN, NA, NE scattering data better than the earlier 
QDCSM calculations, which fit the data only qualitatively [JTI, I2| . Consequently, we expect 
this extended QDCSM calculation will provide a better model estimate of the mass of the 
di-tt. 
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In the following section, we give a short description of the extended QDCSM. The results 
and a discussion are presented in section III. 



II. SHORT DESCRIPTION OF THE EXTENDED QDCSM 



The details of the QDCSM and its extension can be found in Refs. [IT[ \L2\ 18 1 and the 



resonating-group calculation method (RGM) has been presented in Refs. |20|, Here we 
present only the (complete) extended model Hamiltonian, wave functions and the necessary 
equations used in the current calculation. 

The Hamiltonian for the 3-quark system is the same as the usual quark potential model. 
For the six-quark system, it is assumed to be 

H, = E(m 1 + i)-T C M+ E ( W«) + Vo(r y ) + V x (r«)) , 
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if i,j occur in different baryon orbits, 
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where XI and \{ are, respectively, the spin, color and flavor operators of the i th quark. 
The spin operators are represented by the three 2x2 Pauli matrices, and the color and 
flavor operators are represented by the eight 3x3 Gell-Mann matrices. The /x 7 are the 
masses of the Goldstone bosons (7 = ir,K,r)). We assume a single quark-meson coupling 
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constant gf/47r for all mesons. The confinement potential V con f{rij) has been discussed in 



Refs.flJ, |2C 

We obtained the values of b, a s and a c by reproducing the N — A mass difference, the 
nucleon mass and applying a stability condition. The mass of the u, d quark is assumed 
to be 1/3 of the nucleon mass to meet the requirements of the resonating group method 
(RGM) calculation while the strange quark mass m s is determined by an overall fit to the 
baryon octet and decuplet. A flavor symmetric overall octet quark- meson coupling constant 
is assumed for all of the octet mesons (it, K, 77). 

The color screening parameter /1 is determined by the binding energy of the deuteron. 
The full set of parameters reproduces other properties of deuteron as well: the size and 
the D-S wave mixing amplitude. It is worth emphasizing that the extended QDCSM has 
only one parameter, the color screening constant /1, that needs to be adjusted to obtain 
the correct binding energy of the deuteron, and that this is sufficient to reproduce well 
all of the properties of the deuteron. The tensor force of pion exchange plays a vital role 
in reproducing the deuteron properties quantitatively, but the results are not extremely 
sensitive to the meson exchange cutoff parameter, tq. 

After introducing Gaussian functions with different reference centers Si, i=l...n, which 
play the role of the generating coordinates in this formalism, and including the wave function 



for the center-of-mass motionPS, the ansatz for the two-cluster wave function used in the 
RGM can be written as 

n „ JQ 3 6 

^ = ^EE E -£lli, a &,e)IlM-3i,e) 

k i=l L fc =0,2 V 4vr a=1 p =A 

[[^ lfcSlfe ( J B lfe )r ?/2fc52fc ( J B 2fc )]^F^(4)] J [x c (5i)x c ( J B2)] W , (2) 



where k is the channel index. For the di-f2, we have k = 1,2, corresponding to the channels 
flfl S=0, L=0 and S=2, L=2. We consider only these channels as they are the only ones 
coupled by the tensor interaction. 

The delocalized single-particle wave functions used in the QDCSM are 



M&,e) = (0 a (^) + e0 a (-^))/iV(e), 



N(e) = v / l + e 2 + 2ee^ 2 / 4fe2 . (3) 
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Ml) = (i) e-*"-^ 

The delocalization parameter, e, is determined by the dynamics of the quark system rather 
than being one of the fixed (adjusted) parameters. 

From the variational principle, after variation with respect to the relative motion wave- 
function x(R) > one obtains the RGM equation 

J H(R, R') X (R')dR' = E J N(R, R')x(R')dR', (4) 

With the above ansatz, the RGM equation, Eq.([|), becomes an algebraic eigenvalue 
equation, 

E c^ L Xf k,ALk = E H c jAL X'f k ' (5) 

j,k,L h j 

k' V k Lfc k' L' 

where N { J k \H i ^ ' ' k ' are the (Eq.(|2])) wave function overlaps and Hamiltonian matrix 
elements, respectively. By solving the generalized eigen value problem, we obtain the energies 
of the 6-quark system and their corresponding wave functions. 



III. RESULTS AND DISCUSSION 



Our model parameters, which have been fixed by matching baryon and deuteron proper- 
ties and using only tt meson exchange, are given in Table I. For comparison, the results of 



our earlier calculation [18j without ir meson exchange are also included in Table I. In order to 
study the dependence of the result on the short range cutoff radius, we choose three typical 
values: 0.6, 0.8 and 1.0 fm, for the short-range cutoff of the ir meson exchange potential. 
This cutoff is necessary in our model approach because all short and intermediate range 
effects are already represented by the quark delocalization and color screening: The short 
range repulsion of the NN interaction is provided by a combination of the color magnetic 
interaction due to gluon exchange and the Pauli principle enforced by the quark structure 
of the nucleon; the intermediate range attraction conventionally modeled by heavy meson 
and multimeson exchange is also provided by quark delocalization and color screening. If 
the contact term of single 7r meson exchange, 8(f), or its smeared version, exp(— Ar)/(Ar), 



were used, double counting[12] would occur. In all cases, we found that the contribution 
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of 7r meson exchange to the dibaryon mass is not large (about 10 MeV or less), primarily 
due to the effect of the short-range cutoff; correspondingly, the deviations from the original 
QDCSM parameters are also quite small. 

Table I. Model parameters and results calculated for di-Q. 





includin 


g 7i meson exchange 


excluding meson exchange 




tq = 0.6 fm 


tq = 0.8 fm 


ro = 1.0 fm 


dynamical cal. 


adiabatic cal. 


m d ,m s (MeV) 


313 , 634 


313 , 634 


313 , 634 


313, 634 


313 , 634 


b (fm) 


0.6010 


0.6015 


0.6021 


0.6034 


0.6034 


a c (MeV fm" 2 ) 


25.40 


25.14 


25.02 


25.13 


25.13 


a s 


1.573 


1.5585 


1.550 


1.543 


1.543 


/j, (fm -2 ) 


0.75 


0.85 


0.95 


1.0 


1.0 


Mass sc (MeV) 


3290.3 


3298.2 


3303.5 


3312 


3350 


Mass cc (MeV) 


3290.2 


3297.9 


3303.2 







sc refers to a k = 1 single channel calculation, cc refers to the coupled channel calculation. 

Comparing the last two columns in Table I, where the n meson exchange tails have not 
been included, we see that the dynamic calculation modifies the results of the previous 
adiabatic calculation significantly. The mass of the flfl, which is close to the threshold 
(3344.9 MeV) in the adiabatic approximation, is reduced below the threshold by about 
33 MeV. Several effects contribute to this difference, but the most important one is that 
the resonating group method calculates the relative motion between the two quark clusters 
rigorously, greatly improving over the rough estimate of the kinetic energy by the zero-point 
oscillation energy used in the adiabatic approximation. 

We note that one should not expect all baryon-baryon (B-B) states to be modified signif- 
icantly by the dynamical calculation. If the effective B-B interaction has a narrow minimum 
with respect to the separation variable, then the zero-point oscillation energy provides a 
good approximation to the kinetic energy of the relative motion of the two clusters. For 
example, the d* mass does not change very much between the adiabatic and dynamic cal- 
culations. However, if the effective B-B interaction is flat, then the zero-point oscillation 
energy may be not a good approximation for the kinetic energy of the relative motion; the 



di-f2 is such a case|L3]. 

The first three columns show the dynamical mass of the di-VL with (only) 7r meson ex- 
change included but with different cutoff parameters. Comparing these three results with 



the fourth column, one can see that the di-f2 mass is reduced by between 9 and 22 MeV. 
This is due to the fact that the parameters of the model are modified by inclusion of the 
pion even though the pion is not exchanged between s quarks. The binding energies vary 
from 55 MeV to 40 MeV with different cutoffs. The larger the cut off, the smaller is the 
effect on the di-fl mass. 

The last two rows show that the effect of channel coupling caused by the meson tensor 
force is quite small here. The S and D wave coupled channel masses (cc) are almost the 
same as those of the single channel approximation (sc). In the deuteron, the effect on the 
mass of the state, although almost an order of magnitude larger and more significant there, 
was also found to be small. By contrast, as shown below, the D-S wave mixing in both cases 
is comparable. 

The quark spin-orbit interaction has been neglected in the Hamiltonian (1). We would 
expect it to have a similarly significant effect on the D-wave channel. However, its effect 
on the di-Q mass is expected to be small because the D-wave channel coupling itself only 
affects the di-fi mass slightly. 

Fig.l shows the relative motion wave functions for two different cutoff parameters. They 
are not sensitive to the short-range cutoff parameter, r . Fig. 2 shows the effective Q — Q 
potential corresponding to the cutoff r = 0.8/m. The other cutoffs give almost indistin- 
guishable effective potentials. 

As has been mentioned above, the effects of short and intermediate range interaction 
arising from 7r meson exchange are represented by quark derealization and color screening, 
so the only long range effect in our model comes from ir exchange. The effects of heavy 
pseudoscalar meson (K, rj) exchange on the mass of the di-f2 have been checked with a 
representative cutoff of r = 0.8/m. The results are presented in Table II. 
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Table II. Mass of di-Q for tq = 0.8/m and full octet pseudoscalar meson exchange. 





only 7r-exchange 


7r,K,T] -exchange 


m u ,m s 


313,634 


313,634 


b(fm) 


0.6015 


0.6022 


a c (MeV fm" 2 ) 


25.14 


25.03 


a s (fm- 2 ) 


1.5585 


1.555 


n 


0.85 


0.9 


Mass Deuteron (MeV) 


1876 


1875.8 




2.0 


1.92 




4.53% 


4.92% 


Mass Di _ n (MeV) 


3298.2 


3300.0 



Repeating the same process mentioned above to fix the parameters, we found that adding 
K and ^-exchange affects the parameters (b, a c , a s ) only slightly; the properties of the 
deuteron can be reproduced just as well by a very small readjustment the color screening 
parameter, /x. Comparing the results in Table I and II, we find the mass of the di-Q is almost 
unaffected by the addition of K and rj exchange. This confirms our expectation that heavier 
(shorter ranged) meson exchanges are already represented by quark effects in our model, so 
that explicit representation of such exchanges beyond a cut-off scale are not important in 
our approach. 

In summary, within the framework of the extended QDCSM, using a resonating group 
method, we estimate the mass of the di-Q to be about 45 MeV lower than the Q — Q 
threshold. This mass is almost unaffected by channel coupling, the value of the short range 
cutoff length, r , and heavier pseudoscalar meson exchanges. The inclusion of single pion 
or heavier meson exchange has only a minor effect on the di-Q mass. These results are 
consistent with the calculations of Kopeliovich and of Zhang et al.El, [23 . 



We should, however, take note that the precision of any prediction of the di-f2 mass is 
limited by that of the model itself. For example, in our model approach, the mass of the Q 
itself is found to be 1651 MeV, or 21 Mev less than the measured value. It could well be 
argued that this theoretical mass of the Q should be used to calculate the threshold. We 
would then conclude that the di-Q mass is no more than 5 MeV below threshold. From 
the model deviations from experimental values for single baryon masses, we estimate a total 
systematic uncertainty of order 30 MeV as a reasonable value to assign to our model. 



Moreover, it is difficult to estimate how much of the model uncertainty is carried over to 
the dibaryon calculation for such a high strangeness quark system, because the model param- 
eters are mainly determined by the experimental data of the nonstrange sector. Therefore, 
when designing a detector system to hunt for the di-fl state, it is theoretically advisable to 
take into account the possibility that the di-fl may be either stable or unstable with respect 
to the strong interaction. 

This work is partly supported by the HIRFL-CSR theoretical center, the NSF and SED 
of China, and partly by the US Department of Energy under contract W-7405-ENG-36. 
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Figure Captions 

Fig. 1. Cluster relative motion wave function of di-Q. 
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Fig. 2. S-wave f2 — f2 effective potential. 
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FIG. 1: Cluster relative motion wave function of di-O. 
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FIG. 2: S-wave £1 — Q effective potential. 
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